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ABSTRACT 

Context. The super-Eddington wind scenario has been proposed as an alternative way for pro¬ 
ducing type la supernovae (SNe la). The super-Eddington wind can naturally prevent the carbon- 
oxygen white dwarfs (CO WDs) with high mass-accretion rates from becoming red-giant-like 
stars. Furthermore, it works in low-metallicity environments, which may explain SNe la observed 
at high redshifts. 

Aims. In this article, we systematically investigated the most prominent single-degenerate 
WD-tMS channel based on the super-Eddington wind scenario. 

Methods. We combined the Eggleton stellar evolution code with a rapid binary population syn¬ 
thesis (BPS) approach to predict SN la birthrates for the WD+MS channel by adopting the super- 
Eddington wind scenario and detailed mass-accumulation efficiencies of H-shell flashes on the 
WDs. 

Results. Our BPS calculations found that the estimated SN la birthrates for the WD-t-MS chan¬ 
nel are ~0.009-0.315x10“^ yr“* if we adopt the Eddington accretion rate as the critical accretion 
rate, which are much lower than that of the observations (<10% of the observed SN la birthrates). 

This indicates that the WD+MS channel only contributes a small proportion of all SNe la. The 
birthrates in this simulation are lower than previous studies, the main reason of which is that new 
mass-accumulation efficiencies of H-shell flashes are adopted. We also found that the critical 
mass-accretion rate has a significant influence on the birthrates of SNe la. Meanwhile, the results 
of our BPS calculations are sensitive to the values of the common-envelope ejection efficiency. 

Key words, binaries: close - stars: evolution - supernovae: general 

1. Introduction 

Type la supernovae (SNe la) have been successfully used as standard candles for measuring cos¬ 
mological distances, resulting in the determination of the accelerating expansion of the Universe 
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that is driven by dark energy (e.g., Riess et al. 1998; Perlmutter et al. 1999). There is a theo¬ 
retical consensus that SNe la are thermonuclear explosions of carbon-oxygen white dwarfs (CO 
WDs) (see Hoyle & Fowler 1960; Nomoto et al. 1997). However, SN la progenitor systems and the 
physics of their explosion mechanism are still uncertain (for recent reviews see Wang & Han 2012; 
Hillebrandt et al. 2013; Maoz et al. 2014; Ruiz-Lapuente 2014), which may affect the precision of 
results of the current cosmological model (e.g., Howell 2011). 

Depending on the nature of mass donor stars, two ways by which CO WDs in binary systems 
can accrete mass towards the Chandrasekhar limit and then cause SNe la have been proposed, 
i.e., the single-degenerate (SD) model and the double-degenerate (DD) model. In the SD model, 
a CO WD accretes material from a non-degenerate star to increase its mass to the Chandrasekhar 
mass limit, which then leads to an SN la explosion (e.g., Whelan & Iben 1973; Nomoto et al. 
1984). In the DD model, two CO WDs in a close binary are brought together by the gravitational 
wave radiation and merge to result in an SN la explosion (e.g., Webbink 1984; Iben & Tutukov 
1984; Nelemans et al. 2001; Geier et al. 2007; Ruiter et al. 2009; Chen et al. 2012). For recent 
observational constraints to theoretical models see, e.g., Maoz et al. (2014), Wang & Han (2012), 
Liu et al. (2012a) and Wang et al. (2013b). 

The most commonly-considered variant in the SD model is the “WDh-MS channel” (usually 
known as the “supersoft channel”), in which the WD accretes material from a main-sequence (MS) 
star or a slightly evolved sub-giant star (e.g., Hachisu et al. 1996, 1999; Li & van den Heuvel 1997; 
Langer et al. 2000; Han & Podsiadlowski 2004; Wang et al. 2010; Claeys et al. 2014; Ablimit et 
al. 2014). This channel has been theorized to be accompanied by recurrent novae and supersoft 
X-ray sources in the past two decades (see Parthasarathy et al. 2007; van den Heuvel et al. 1992; 
Rappaport et al. 1994). Note that we did not study the SD WD-hRG channel in this paper as the SN 
la birthrate from that channel is far below the birthrate from the WDh-MS channel (e.g., Yungelson 
& Livio 1998; Lii et al. 2006; Mennekens et al. 2010; Claeys et al. 2014). 

The accretion of H-rich material on WDs is crucial for understanding SNe la from the WDh-MS 
channel. If the mass-accretion rate is too high, the WD will evolve into a red-giant-like star due to 
the pileup of the accreted material on its surface. If the mass-accretion rate is too low, the nuclear 
burning is unstable, H-shell flashes occur and expel the accreted shell. Thus, there is only a narrow 
range where the accreted H-rich material could steadily burn on the surface of WDs. This leads to a 
fact that theoretical SN la birthrates from the WDh-MS channel are far below those of observations 
(e.g., Yungelson et al. 1996). The optically thick wind assumption was proposed to solve this 
problem, in which the H-rich material is transformed into He at a critical rate, while the unprocessed 
material is blown off in the form of the optically thick wind (see Hachisu et al. 1996,1999; Nomoto 
et al. 2007). However, the optically thick wind assumption is still quite controversial (see Langer 
et al. 2000). Besides, this assumption does not work when the metallicity is lower than 0.002 (see 
Kobayashi et al. 1998), while, as a matter of fact, SNe la with a metallicity lower than this threshold 
have been observed (e.g., Prieto et al. 2008; Badenes et al. 2009). 

Ma et al. (2013) recently proposed an alternative way to replace the optically thick wind 
regimes, which is called the super-Eddington wind scenario. In this scenario, the super-Eddington 
wind is triggered when the luminosity of the accreting WD exceeds the Eddington luminosity. 
The super-Eddington wind can naturally prevent the CO WDs with high mass-accretion rates from 
forming red-giant-like stars. In contrast to the optically thick wind, the super-Eddington wind does 



B. Wang et al.: A super-Eddington wind scenario of SNe la; BPS calculations 3 

not significantly depend on the metallicity, which could be triggered even for Z = 10 ® (see Fig. 4 
of Ma et al. 2013). This indicates that this scenario may have a contribution to SNe la at high red- 
shifts. In this work, we aim to investigate theoretical SN la birthrates and delay-time distributions 
of the WDh-MS channel by adopting the super-Eddington wind scenario. 

As mentioned above, at low mass-accretion rates, the accreting WD would experience H-shell 
flashes like nova outbursts. However, the H-shell flashes have not been well understood due to 
the numerical difficulties in nova outbursts. For example, only one cycle of H-shell flash has been 
followed in most simulations of nova outbursts; multicycle evolution of the H-shell flash is difficult 
to compute (e.g., Kovetz & Prialnik 1985; Iben et al. 1992; Schwartzman et al. 1994). The H-shell 
flashes are important for understanding the SN la progenitor scenario studied in this paper; a low 
mass-accumulation efficiency of H-shell flashes will lead to a low birthrate of SNe la. Denissenkov 
et al. (2013a) recently constructed multi-cycle nova evolutionary sequences with WDs using the 
stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA). Adopting a similar 
method to Denissenkov et al. (2013a), Ma et al. (2013) obtained the mass-accumulation efficiency 
of H-shell flashes for a wide range of WD masses. A more detailed binary population synthesis 
(BPS) approach is needed to study the influence of the H-shell flashes on the hnal results and to 
give SN la birthrates, which is also addressed in this work. 

This paper is organized as follows. The numerical code for the binary evolution calculations 
is described in Sect. 2, and the corresponding results are provided in Sect. 3. The BPS method and 
results are presented in Sect. 4. Our discussion and summary are then presented in Sect. 5. 


2. Numerical code for binary evolution 


In the WD-hMS channel, the mass donor star is a MS star or a slightly evolved sub-giant star, which 
transfers H-rich material on the surface of the WD. The accreted H-rich material is burned into He, 
and then the He is converted into carbon and oxygen, resulting in the mass increase of the WD. We 
assume that the WD explodes as an SN la once its mass grows to 1.378 M©, which is the critical 
mass limit of non-rotating WDs for carbon ignition (e.g., Nomoto et al. 1984)Q We calculated the 
evolution of the WD-i-MS systems using the Eggleton stellar evolution code (Eggleton 1973; Han 
et al. 1994; Pols et al. 1995, 1998; Eggleton & Kiseleva-Eggleton 2002). The initial setup and basic 
input physics for this code are similar to those in Wang et al. (2010). The MS star models in our 
calculations are composed of metallicity Z = 0.02, H abundance X - 0.70, and He abundance 
Y = 0.28. 

The mass transfer begins once the MS star hlls its Roche-lobe. We dehne the values of mass- 
accumulation efficiency for H-shell burning on the surface of the WD, t/h, as 


MEdd/lMll, IMll > MEdd, 


7h = i 1, 


MEdd ^ \M2\ > Mst, 


IM2I < 


( 1 ) 


where M 2 , MEdd and Mjt are the mass-transfer rate, the Eddington critical accretion rate, and the 
minimum accretion rate for stable H-shell burning, respectively. In the following, we make a de- 

' We did not consider the influence of rotation on the evolution of the accreting WD in this work; the 
maximum stable mass for a rotating WD may be above the standard Chandrasekhar mass limit (see Uenishi 
et al. 2003; Yoon et al. 2004; Wang et al. 2014). 
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Fig. 1. Properties of H-shell burning on the surface of accreting WDs in the plane of WD mass and 
mass-accretion rate. A comparison of the critical mass-accretion rate between the super-Eddington 
wind scenario and the optically thick wind model was made. The solid lines are from Ma et al. 
(2013), while the dashed lines are from Nomoto et al. (2007). 

tailed description for the above equation: (1) If the mass-transfer rate, M 2 , is above the Eddington 
critical accretion rate, Mndd, we assume that the accreted H steadily burns on the surface of the 
WD, and that the H-rich material is converted into He at a rate of 

Medd = 5.975 x 10’®(M^p - 3.496 M^d -h 4.373M^p - 2.226 Mwd + 0.406), (2) 

where Mwd is the mass of the WD in units of M©, and Mndd is in units of Moyr^'. The unprocessed 
material can be lost from the system in the form of the super-Eddington wind with a mass-loss rate 
of (IM 2 I - MEdd)- In Eig. 1, we make a comparison of the critical mass-accretion rate between the 
super-Eddington wind scenario and the optically thick wind model (for more discussions about this 
hgure see Sect. 5). (2) When IM 2 I is smaller than Mndd but higher than the minimum accretion rate 
for stable H-shell burning, the burning of the H layer is steady and no mass is lost from the 
system. The values of Mjt are 

Mst = 3.057 X 10''^(M^o - 0.386Mwd + 0.027), (3) 

where Mst is in units of Meyr"^ (see Eig. 1). (3) H-shell flashes are triggered once IM 2 I is lower 
than Mst- The mass-accumulation efficiency for H-shell flashes, 77 ^, is linearly interpolated from 
a detailed grid (see Eig. 2), in which a wide range of WD masses and mass-accretion rate were 
calculated by Ma et al. (2013). 

Helium is assumed to be ignited once the mass of the He layer under the H-shell burning reaches 
a certain value. The mass-growth rate of the He layer is defined as 

Mho = (4) 

If He-shell flashes occur, we assume that a part of the layer mass is blown off. In this case, we adorn 
the prescription of Kato & Hachisu (2004) for mass-accumulation efficiency of He-shell flasheso 
We define the mass-growth rate of the WD, Mco. as 

Mco - ?7HeAfHe = ?7He?7H|Af2L (5) 

^ Yoon et al. (2004) argued that the burning of He layer is less violent when rotation is considered, which 
may increase the mass-accumulation efficiency of He-rich material. 
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Fig. 2. The mass-accumulation efficiency of H-shell flashes, t^h, is plotted against the mass- 
accretion rate for various initial WD masses. The data points are taken from Ma et al. (2013). 




Fig. 3. A representative example of binary evolution calculations, where the binary system is in the 
weak H-shell flash phase at the moment of the SN explosion. In panel (a), the solid, dashed and 
dash-dotted curves present M 2 , Mco and Mwd, varying with time, respectively. In panel (b), the 
evolutionary track of the donor in the Hertzsprung-Russell diagram is shown as a solid curve, and 
the evolution of orbital period is shown as a dash-dotted curve. Dotted vertical lines in both panels 
indicate the position where the WD produces an SN la. 


where //He is the mass-accumulation efficiency of He-shell flashes. 

The evolution of the orbital periods of these WD-i-MS systems is mainly decided by the wind 
mass-loss of the WD, which carries away the orbital angular momentum of the binaries. We assume 
that the mass lost from these systems takes away the specific orbital angular momentum of the WD, 
but the mass lost from the wind of the mass donors is negligible (e.g., Wang et al. 2010). Finally, 
dense model grids of binary evolution calculations are obtained. 


3. Binary evolution results 

In Fig. 3, we present a representative example of binary evolution calculations. The binary system is 
in the weak H-shell flash phase at the moment of SN explosion. The initial condition for this binary 
is (Mj, log(PVday)) = (2.1, 0.9, 0.4), where M^, M\^^ and P' are the initial masses of the 


log P (days) 
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Fig. 4. Regions in the logP’ - plane for WD-i-MS binaries that ultimately produce SNe la for 
various initial WD masses, where and M'^ are the initial orbital period and the initial mass of the 
MS star, respectively. 


MS star and the CO WD in solar masses, and the initial orbital period in days, respectively. The MS 
star firstly fills its Roche-lobe after the exhaustion of central H (it now contains a He core), which 
leads to early Case B mass transfer; the orbital period of the binary system is not changed until the 
MS star hlls its Roche-lobe, see panel (b) of Fig. 3. The mass-transferrate IM 2 I exceeds Mndd soon 
after the onset of Roche-lobe overflow, resulting in a wind phase in which a part of the transferred 
mass is blown off in the form of the super-Eddington wind, and the remainder is accumulated on 
the WD; the wind carries away the specific orbital angular momentum of the WD, which results 
in the decrease of the orbital period. After about 1.6 x 10^ yrs, IM 2 I drops below Meh, but is still 
higher than Mjt- At this stage, the H-shell burning is stable, but the super-Eddington wind stops. 
With the decrease of IM 2 I. the WD quickly enters a weak H-shell flash phase. The orbital period 
of the binary starts to increase once the accreting WD is more massive than the donor star. The 
WD always grows in mass until it explodes as an SN la in the weak H-shell flash phase. When the 
WD increases its mass to 1.378 Mq, the mass of the donor star is = 1.1168 M© and the orbital 
period is log(f’™/day) = 0.2092. 

Eigure 4 presents the initial contours for producing SNe la in the orbital period and secondary 
mass plane for various initial WD masses (i.e., - 0.90, 1.0, 1.1 and 1.2 M©). The enclosed 

region almost vanishes for WDs with 0.9 M©, which is then assumed to be the minimum initial WD 
mass for producing SNe la from the WDh-MS channel. We note that the critical mass-accretion rate 
for WDs is still uncertain (e.g., Cassisi et al. 1998; Danger et al. 2000; Nomoto et al. 2007; Shen & 
Bildsten 2007; Bours et al. 2013). In the observations, some binaries which contain a neutron star or 
a black hole can lose their masses at rates that exceed the Eddington accretion rate of the compact 
object by several times (e.g., Cyg X-2 and SS 433; Podsiadlowski & Rappaport 2000; Blundell et 
al. 2001). In Eig. 5, we show the initial contours for producing SNe la after adopting 2MEdd as the 
critical mass-accretion rate; we did not change the values of the minimum accretion rate for stable 
H-shell burning (i.e., Mjt) in our calculations. This is an unphysical extreme assumption, but we 
can use this to examine the influence of the critical mass-accretion rate on the final results. Erom 







B. Wang et al.: A super-Eddington wind scenario of SNe la; BPS calculations 


7 



log (days) 

Fig. 5. Similar to Fig. 4, but for results obtained after adopting 2MEdd as the critical mass-accretion 
rate. 

this figure, we can see that the initial contours for producing SNe la are enlarged significantly. This 
is because a larger regime for stable H-shell burning is adopted. 

The left boundaries of the initial contours in Figs 4-5 are constrained by the condition that 
Roche-lobe overflow begins when the mass donor star is in the zero-age MS stage, while sys¬ 
tems beyond the right boundary experience mass transfer at a relatively high rate due to the rapid 
expansion of the donor stars in the Hertzsprung gap and they lose too much mass via the super- 
Eddington wind, preventing the WDs from increasing their masses to the Chandrasekhar mass. 
The upper boundaries are set by a high mass-transfer rate owing to a large mass ratio between the 
mass donor star and the WD, while the lower boundaries are constrained by the condition that the 
mass-transfer rate should be high enough to ensure the mass growth of the WD during H-shell 
flashes. 

4. BPS assumptions and results 

4.1. Method of BPS 

In order to obtain SN la birthrates and delay times for the WD-hMS channel, a series of Monte 
Carlo BPS simulations are performed. The random number generator algorithm in our Monte Carlo 
simulations is from Press et al. (1992), which provides a uniform probability density in the range 
of (0, 1). The repetition period of the random number generator in our simulations is >10'^, which 
is infinite for practical calculations. 

For each BPS simulation, we employed the Hurley rapid binary evolution code (Hurley et al. 
2002) to simulate the evolution of 10^ primordial binaries, and adopted a metallicity of Z = 0.02. 
These binaries are followed from the star formation to the formation of the WD-hMS systems based 
on three binary evolutionary ways (see Fig. 1 of Wang & Han 2012). If the initial parameters of a 
CO WD-hMS system at the onset of the Roche-lobe overflow are located in the SN la production 
regions in the plane of (logP', for its specific then an SN la is assumed to occur. We 

adopt linear interpolation if is not among the masses listed in Figs 4-5. 

Seven sets of Monte Carlo BPS simulations are conducted to investigate the birthrates of SNe la 
(see Table 1), in which sets 1 and 2 are two best models from Camacho et al. (2014) that result in a 
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Table 1. Galactic SN la birthrates for seven simulation sets, 
in which sets 1 and 2 are two best models from Camacho et 
al. (2014). Notes: acE - CE ejection efficiency; n(q) = initial 
mass-ratio distribution; IMF = initial mass function; - 
critical mass-accretion rate; v = Galactic SN la birthrate. 


Set 

UCE 

n(q) 

IMF 

Mcr 

V (10 ^ yr ') 

1 

0.3 

1 

MS79 

Affidd 

0.315 

2 

0.3 

OC q 

MS79 


0.137 

3 

1.0 

1 

MS79 

Affidd 

0.110 

4 

1.0 

1 

S86 

Mfidd 

0.070 

5 

1.0 

1 

MS79 

2MEdd 

0.955 

6 

3.0 

1 

MS79 

3fEdd 

0.009 

7 

3.0 

1 

MS79 

2MEdd 

0.712 


good overall fit to the properties of the WD-hMS binary population in the Sloan Digital Sky Survey 
(SDSS) Data Release 7. The following basic assumptions for these Monte Carlo simulations are 
adopted: 

(1) A constant star formation rate (SFR) of 5 Moyr^' is simply assumed over the past 14 Gyrs 
or, alternatively, it is modeled as a delta function, that is a single instantaneous starburst (10*® M© in 
stars is assumed). We intend the constant SFR to approximate spiral galaxies, and the delta function 
to provide a rough description of elliptical galaxies or globular clusters (e.g., Wang et al. 2013a). 

(2) All stars are assumed to be members of binaries. The primordial binary samples are gen¬ 
erated via a Monte Carlo method and an initially circular orbit was supposed for all binaries. The 
distribution of initial orbital separations is assumed to be constant in log a for wide binaries, where 
a is the orbital separation (e.g., Han et al. 1995). 

(3) The initial mass function (IMF) of the primordial primary star is taken from Miller & Scalo 
(1979, MS79). Alternatively, we also consider the IMF of Scalo (1986, S86). 

(4) A constant mass-ratio distribution is taken, i.e., n{q) - 1, or a distribution proportional to 
the mass ratio, n(q) oc q (e.g., Camacho et al. 2014). 

(5) The common-envelope (CE) evolution is still uncertain (e.g., Zuo &. Fi 2014). The standard 
energy equations are used to describe the output during the CE stage (e.g., Webbink 1984). There 
are two uncertain parameters for this prescription of the CE ejection, i.e., acE (the CE ejection 
efficiency) and A (a structure parameter that depends on the evolutionary stage of the mass donor). 
A is usually set to be 0.5 for the purpose of constraining acE (see de Kool 1990). Here, we set acE 
to be 0.3, 1.0 and 3.0 to examine its influence on the final results. 

4.2. Results of BPS 

Based on the seven sets of simulations for the WD-hMS channel, the estimated birthrates of SNe la 
are strongly sensitive to the choice of some initial conditions, e.g., the critical mass-accretion rate, 
the CE ejection efficiency, initial mass-ratio distribution and initial mass function, etc. Significantly, 
if we adopt 2MEdd as the critical mass-accretion rate, the SN la birthrate will increase to be a 
maximum value. This is because a larger initial contours for producing SNe la are adopted. 
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Fig. 6. Evolution of SN la birthrates for a constant star-formation rate (SFR = S.OMoyr ') with 
different BPS simulation sets. 



log(i/yr) 


Fig. 7. Delay-time distributions of SNe la with different BPS simulation sets. The open circles are 
taken from Maoz et al. (2011), while the filled squares are from Totani et al. (2008). 


In Fig. 6, we show the evolution of SN la birthrates with time by adopting metallicity Z = 0.02 
and star-formation rate SFR = 5Moyr“’. This study presents an SN la birthrate of ~0.009-0.315 x 
10“^ yr“' by adopting MEdd as the critical accretion rate, which is far below the observed value (i.e., 
3 -4 X 10“^ yr Cappellaro & Turatto 1997)0Thus, we speculate that the WDh-MS channel only 
contributes to a small part of all SNe la (<10%); some other formation channels or mechanisms 
may also have a contribution to SNe la (for a recent review see Wang & Han 2012). If we adopt 
2MEdd as the critical accretion rate, the birthrate of SNe la will increase to ~0.712-0.955 x 10“^ yr“' 
(see sets 5 and 7). This value is still lower than that of observations, but it provides an upper limit 
on the birthrate of SNe la for the WD-i-MS channel. This means that the variation of the critical 
accretion rate has a significant influence on the SN la birthrates. 

SN la delay times are defined as the time interval from star formation to the SN explosion. 
In Fig. 7, we display the SN la delay-time distributions for the WD-hMS channel with different 

^ The upper mass limit for CO WDs has been suggested to be ~1.07 Mq (see Umeda et al. 1999). If this 
mass limit is adopted in Fig. 4, the SN la birthrate from the WD+MS channel will decrease to ~0.009-0.215 x 
10-3 yr-‘. 
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Fig. 8 . Similar to Fig. 7, but for the case of set 1, where ten random initial seeds are adopted in our 
simulations. 




Fig. 9. Distribution of properties of initial WD-hMS systems that can ultimately produce SNe la. 
Panel (a) shows the results in the logP’ - plane, while panel (b) gives the results in the 
log P‘ - plane, where P\ and are the initial orbital period, the initial WD masses, and 
the initial mass of the MS star, respectively. Here, we present the case of set 1. 


BPS simulation sets, which are obtained from a single starburst. From this figure, we can see that 
the estimated delay times for this channel are ~110Myrs-1200Myrs after the starburst, which 
may have a contribution to the SNe la with intermediate delay times. Note that a low value of 
acE results in a higher birthrate of SNe la. This is because, for a low value of acE, the primordial 
binary systems need to release more orbital energy to eject the CE, and then produce WD-hMS 
binaries. In this case, WD-hMS binaries tend to have slightly closer orbital separations, and are 
more likely to be located in the SN la production region (e.g.. Figs 4-5). We note that, in any 
Monte Carlo simulations the initial seed of the random number generator results in a given and 
specific simulated population. In Fig. 8, we show the delay-time distributions of SNe la but for the 
case of set 1, in which ten random initial seeds are adopted in our simulations. From this figure, 
we found that there is no significant difference among these simulations, i.e., our Monte Carlo BPS 
results are not sensitive to the initial seed of the random number generator. 

In the observations, some WD-tMS systems are suggested to be candidates of SN la progenitors 
(for recent reviews see Wang & Han 2012; Parthasarathy et al. 2007 and references therein). In 
Fig. 9, we give some properties of initial WDh-MS systems that can produce SNe la based on 
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our BPS calculations, which would be useful for identifying potential progenitor candidates of 
SNe la. For these WDh-MS systems, SN la explosions in our simulations occur for the ranges 
M^d ~ 0.9-1.2Mo, ~ 2.1-3.1 Mq and P* ~ 0.5-5days. However, recent observations in the 

SDSS have not found massive MS stars in WDh-MS systems, where the MS stars have their masses 
<0.6 Mq (e.g., Zorotovic et al. 2011; Nebot Gomez-Moran et al. 2011; Rebassa-Mansergas et al. 
2012, 2013; Toonen & Nelemans 2013; Camacho et al. 2014; Li et al. 2014), which cannot afford 
enough material for the WDs to increase their masses to the Chandrasekhar limit. 


5. Discussion and conciusions 

Previous studies usually assumed that, when the mass-transfer rate is lower than jMcr but higher 
than |Mcr (see dashed lines in Fig. 1), a weak H-shell flash is triggered but no mass is lost from 
the binary system, i.e., the H is totally converted into He during the mass-accretion phase (e.g., 
Hachisu et al. 1999; Han & Podsiadlowski 2004; Meng et al. 2009; Chen & Li 2009; Wang et 
al. 2010; Tauris et al. 2013). However, in this paper we adopted a detailed mass-accumulation 
efficiency for H-shell flashes once the mass-transfer rate is below Mjt (see Fig. 2), the value of 
which is lower than previous studies, leading to a low birthrate of SNe la. Note that there is a 
heated debate about the mass-accumulation efficiency on WDs; it is still uncertain whether nova 
outbursts remove more mass from the WD than the accreted material (e.g., Prialnik & Kovetz 1995; 
Cassisi et al. 1998; Yaron et al. 2005; Idan et al. 2012; Newsham et al. 2013). 

Denissenkov et al. (2013b) recently proposed a new type of WDs, i.e., hybrid CONe WDs, 
which have an unburned CO core surrounded by a thick ONe shell. In a subsequent investigation, 
Chen et al. (2014) found that this type of hybrid WDs could be as large as 1.3 M©. According 
to the works of Chen et al. (2014), Meng & Podsiadlowski (2014) studied the hybrid WDh-MS 
channel of SN la progenitors based on a detailed BPS approach, and gave the birthrates ~0.026 - 
0.24 X 10“^ yr“^. If this kind of hybrid WDs can really undergo a thermonuclear runaway when 
increasing their masses to a critical limit, then the SN la birthrates from the WDh-MS channel can 
be increased significantly. Note that hybrid WDs could reach a state of explosive carbon ignition, 
but this depends on some mixing assumptions and the convective Urea process (see Denissenkov 
et al. 2015). 

When the WD in the WDh-MS channel explodes as an SN la, the mass donor star in this channel 
would survive and potentially be identifiable (e.g., Liu et al. 2012b; Pan et al. 2014). Tycho G was 
referred as the surviving companion of Tycho’s SN from the WDh-MS channel (e.g., Ruiz-Lapuente 
et al. 2004; Wang & Han 2010), although it is still debatable (e.g., Ihara et al. 2007; Kerzendorf et 
al. 2009; Liu et al. 2013). The surviving companions from this channel may also provide a way to 
explain the formation of single low-mass He WDs in the observations (e.g., Justham et al. 2009; 
Wang & Han 2010). 

The process of mass accretion on CO WDs is important for understanding the SD model of SNe 
la, which has a direct effect on the theoretical birthrates of SNe la (see alsoBours et al. 2013). In this 
article, by employing the Eggleton stellar evolution code and adopting detailed mass-accumulation 
efficiencies of H-shell flashes on the WDs, we performed detailed binary evolution calculations for 
the WDh-MS channel based on the super-Eddington wind scenario, and gave the initial parameter 
spaces for producing SNe la with various initial WD masses. We found that a CO WD with its mass 
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as low as 0.9 Mq can potentially accrete enough mass and then reach the Chandrasekhar mass limit 
if we adopt Medd as the critical accretion rate. By incorporating the results of the detailed binary 
evolution calculations into a detailed BPS approach, we obtained the evolution of SN la birthrates 
with time. The SN la birthrate from the WD-i-MS channel is much lower than the value obtained 
from observations. This means that the contributions of the WD-hMS channel to the total SN la 
birthrates are quite inefficient. However, this result is strongly sensitive to uncertainties in some 
initial conditions, especially the efficiency of mass accretion. If we adopt 2MEdd as the critical 
accretion rate, the SN la birthrate is increased significantly. Additionally, the values of acE also 
have an influence on the birthrates of SNe la. To set further constraints on the WD-hMS channel, 
large samples of observed massive WDh-MS systems and surviving companions are needed, and 
the process of mass accretion on WDs should be examined carefully. 
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